Superselection from canonical constraints 5 

Michael J. W. Hall 
Theoretical Physics, IAS, 
Australian National University, 
Canberra ACT 0200, Australia 



Abstract 

The evolution of both quantum and classical ensembles may be de- 
scribed via the probability density P on configuration space, its canon- 
ical conjugate S, and an ensemble Hamiltonian H[P, S]. For quantum 
ensembles this evolution is, of course, equivalent to the Schrodinger 
equation for the wavefunction, which is linear. However, quite simple 
constraints on the canonical fields P and S correspond to nonlin- 
ear constraints on the wavefunction. Such constraints act to prevent 
certain superpositions of wavefunctions from being realised, leading 
to superselection-type rules. Examples leading to superselection for 
energy, spin-direction and 'classicality' are given. The canonical for- 
mulation of the equations of motion, in terms of a probability density 
and its conjugate, provides a universal language for describing classical 
and quantum ensembles on both continuous and discrete configuration 
spaces, and is briefly reviewed in an appendix. 
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1 Introduction 



The Schrodinger equation for a quantum system is linear, implying that the 
superposition of any two solutions is also a solution. However, some combina- 
tions of states have never been observed, including coherent superpositions of 
integer and half-integer spins, electric charges, and Schrodinger's cat (wanted 
dead and alive). Possible explanations for why such superpositions are not 
observed fall into two logical categories: 

(i) measurement supers election rules: such superpositions may be allowed, 
but physical limitations on measurement prevent their observation; 

(ii) state supers election rules: such superpositions are not physically allowed. 
Each of these categories may of course be further subdivided into 'in princi- 
ple' and 'in practice' subcategories. 

A well known example of a measurement superselection rule is charge su- 
perselection in quantum field theory. In particular, integration of the opera- 
tor relation V.E = 47rp, relating the electric field and charge density, implies 
that the total charge operator Q = J d 3 x p can be expressed as an integral 
over a surface containing any given finite region, and hence that all quasilocal 
field observables (those vanishing outside such regions) must commute with 
Q [T]. Thus, the assumption that all physical observables are quasilocal im- 
plies that no measurement can distinguish relative phases of superpositions 
of different charge eigenstates, eg, between \qi) + \q 2 ) and \qi) — \q 2 ). Note 
that this is not a state superselection rule, as such superpositions are not in 
themselves forbidden. 

Measurement superselection rules also include environment-induced su- 
perselection, forbidding the observation of coherent superpositions of classical 
apparatus states (and cats), due to decoherence arising from interaction be- 
tween the apparatus and its (in practice unobservable) environment [3 El E] ■ 
It has been suggested that environment-induced superselection rules may also 
provide an alternative basis for explaining charge superselection 

State superselection rules are stronger than (and imply) measurement 
superselection rules (one cannot observe what does not exist). They have 
typically been obtained by appealing to group symmetry arguments. For 
example, assuming invariance of the wavefunction (up to a phase) under 
rotations of 27r, the nontrivial nature of projective representations of the ro- 
tation group implies that superpositions of states of integer and half-integer 
spin are not permitted (HI El IE]- Note that requiring a similar invariance 
for the wavefunctional of a charged field, under global phase shifts of the 
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field, implies that superpositions of states of different total charge are not 
permitted. However, as noted by Weinberg (his italics): u it may or may not 
be possible to prepare physical systems in arbitrary superpositions of states, 
but one cannot settle the question by reference to symmetry principles, be- 
cause whatever one thinks the symmetry group of nature may be, there is 
always another group whose consequences are identical except for the absence 
of superselection rules" jHj (see also Ref. [5]). This emphasises the fact that, 
in the end, superselection rules are incorporated into quantum theory via 
assumptions rather than as strict consequences. 

The aim of this note is to point out an alternative mechanism for state 
superselection rules: nonlinear constraints on the wavefunction. Consider- 
ation of such constraints is perfectly natural - and in principle unavoidable 
- when one places physical significance on the existence of an action princi- 
ple for quantum ensembles. In particular, any Hamiltonian formulation of 
such an action principle requires two canonically conjugate fields on config- 
uration space (eg, ip and ip* [7J [TU], Reip and Im-0 [TT], or \ip\ and arg^ 
[T^l IT!H H3j). Specifying any constraints on these fields is then merely part 
of the general theory of Hamiltonian systems (that underlies, for example, 
classical mechanics, quantization on curved spaces, and standard discussions 
of quantum gauge symmetries [T5]). Note that it is not convincing to argue 
that any such constraints must in fact be linear in ip, on the grounds that the 
set of physical wavefunctions should be closed under superposition, as this 
argument amounts to simply assuming that there are no state superselection 
rules. 

There is another fundamental justification for considering constraints that 
are nonlinear with respect to the wavefunction. In particular, as will be re- 
viewed in section 2, the evolution of both classical and quantum ensembles 
may be described in terms of an action principle for two canonically conjugate 
fields P and S, where P denotes the probability density on the configura- 
tion space. In this common formalism, the primary differences between the 
classical and quantum cases are due to different ensemble Hamiltonians. In 
particular, there is no difference in the treatment of quantum and classical 
constraints: these are expressed in precisely the same way, in terms of the 
canonical fields P and S. Hence, since for classical ensembles one cannot 
even define, much less restrict to, the notion of constraints linear in some 
'wavefunction', it is rather difficult to justify any such a priori restriction for 
quantum ensembles. 

It turns out that quite simple constraints on P and S can lead to inter- 
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esting state superselection rules for quantum ensembles. Moreover, such su- 
perselection rules can be qualitatively rather different from those which have 
been previously considered. For example, they can act to not only restrict 
quantum ensembles to members of a set of orthogonal subspaces in Hilbert 
space, but to further permit only certain types of superpositions within these 
subspaces. Indeed, they can even act to restrict quantum ensembles to mem- 
bers of a set of non-orthogonal states. It should be emphasised that in all 
cases the evolution of the system is still linear, being described by the usual 
Schrodinger equation. 

In the following section the Hamiltonian formulation of the equations of 
motion for classical and quantum ensembles of particles is briefly described. 
Differences between canonical constraints on such ensembles, and Dirac-type 
constraints on the wavefunction, are discussed in section 3. Examples of con- 
straints leading to energy, spin-direction and Gaussian superselection rules 
are given in section 4, and conclusions presented in section 5. 

The Hamiltonian formulation of the description of ensembles on config- 
uration space is of some interest in its own right, being quite universal in 
nature. It is applicable to both continuous and discrete configuration spaces; 
to both classical and quantum ensembles; and even to non-orthogonal repre- 
sentations of quantum ensembles (eg, coherent-state representations). Some 
general features of this formulation are therefore collected in an appendix. 

2 Canonical equations of motion 

The existence of an action principle for the evolution of a statistical ensem- 
ble, in a continuous configuration space, implies equations of motion of the 
Hamiltonian form (see appendix) 

(W_5H_ dS_ _5H_ 

Here P(x, t) denotes the probability density on the configuration space, 
S(x, t) is the quantity canonically conjugate to P, and H[P, S] is the ensem- 
ble Hamiltonian. The notation S/5P, 5/5S denotes functional derivatives 
with respect to P and S respectively (see appendix). 

For example, for a classical ensemble of particles of mass m, acted on by 



4 



an external potential V(x), define the ensemble Hamiltonian 



H C [P,S] := J dxP 



\vs 



2 



+ V(x) 



(2) 



2m 

The equations of motion then follow via Eq. (JTJ) and Eq. (|2~3|l of the appendix 

as 

f + V.(^S/ m) =0, f + ^! + ^0. (3) 

The first of these is a continuity equation, ensuring conservation of proba- 
bility, while the second is the classical Hamilton- Jacobi equation JH]. Note 
that if VS* is interpreted as the momentum of a particle at position x, then 
the classical ensemble Hamiltonian H c corresponds to the average ensemble 
energy. However, such an interpretation will not be relied upon here. 

As a second example, for a quantum ensemble of such particles define the 
ensemble Hamiltonian 

ft 2 r iVlog-Pl 2 
H q [P, S] := H C [P, S} + -j dxP 1 - °* 1 . (4) 

The corresponding equations of motion for P and S follow as 

- + V.(PV5/ m) =0, _ + L_l + v + O = 0, (5) 

where Q := — (k 2 /2m)P~ 1 ^ 2 'V 2 P 1 ^ 2 is the so-called 'quantum potential' [T7| . 
It is well known that if one defines ip := p 1 / 2 e lS / h , these equations of motion 
are equivalent to the Schrodinger equation 

It may be checked that the quantum ensemble Hamiltonian H q is equal to 
the average quantum ensemble energy / dxip*(—h V 2 /2m + V)?p (see also 
appendix). Note also that H q reduces to H c in the classical limit h — > 0. 

These examples indicate that both classical and quantum ensembles may 
be described via two canonically conjugate fields P and S on configuration 
space, with differences in evolution corresponding to differences in the re- 
spective ensemble Hamiltonians H c and H q . Note from Eq. (j3J) that these 
ensemble Hamiltonians differ by a nonclassical kinetic energy term, inversely 
proportional to mass m. Some general properties of the canonical formu- 
lation for ensembles on configuration space (whether classical or quantum, 
continuous or discrete) are discussed in the appendix. 
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3 Dirac constraints vs canonical constraints 



The canonical formulation of the equations of motion provides a common 
basis for discussing both classical and quantum ensembles, and in particular 
for treating classical and quantum constraints on an equal footing. More 
significantly, for the purposes of this note, the canonical formulation in fact 
leads to a more general treatment of quantum constraints than does the 
standard approach introduced by Dirac [TH] . 

In particular, the standard approach relies on replacing classical con- 
straints of the form C(x,p) = by operator constraints on the wavefunction, 
of the form UHl EEE] 

C(x, -V)V> = 0. (7) 

Consistency of these constraints with the equations of motion and with each 
other then leads to a number of secondary constraints of a similar form. What 
is important to note here is that in this approach all constraints are linear 
with respect to the wavefunction. Hence any superposition of two wavefunc- 
tions, each satisfying the constraints, will also satisfy the constraints, i.e., 
the standard approach cannot lead to supers election rules. 

In contrast, in the canonical formulation one is dealing, in both the quan- 
tum and the classical cases, with a Hamiltonian system for two conjugate 
fields, P and S. Hence, any primary constraints are expressed in precisely 
the same manner, 

K[P,S] = 0, (8) 

in both cases ^3]. Any fundamental classical/quantum differences in this 
regard arise because the corresponding ensembles have different ensemble 
Hamiltonians, so that any secondary constraints, following from consistency 
with the equations of motion, will in general be different ^3]. It may be 
assumed that the functional forms of the classical and quantum primary 
constraints are identical, i.e., K q = K c (this is the simplest possible corre- 
spondence). Note, however, that physical consistency only requires identity 
to lowest order in h. 

Now, every Dirac-type constraint on the wavefunction, of the form of 
Eq. (J2J), clearly corresponds to two canonical constraints on the fields P and 
S, of the form of Eq. flSJ) (obtained by substituting ip = p l / 2 e lS / h in Eq. (J7|) 
and taking real and imaginary parts). For example, the Dirac constraint 
(ix V)^> = (rotational symmetry of i/j) maps to the corresponding canonical 
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constraints x x VP = 0, x x VS* = (rotational symmetry of P and S 
respectively). More general symmetries of ip also map to two corresponding 
symmetries for P and S (eg, invariance of ip under rotation by 2rr, up to a 
phase factor, maps to the invariance of P, and to the invariance of S up to 
an additive constant). 

However, the set of possible canonical constraints is much larger than 
those corresponding to Dirac-type constraints, and in particular includes con- 
straints on P and S that correspond to nonlinear constraints on ip (examples 
are given in the next section). Such nonlinear constraints are not invariant 
under superposition, and hence lead directly to state superselection rules for 
quantum ensembles. 

As remarked in the Introduction, the possibility of classical constraints 
which are not linear for the particular combination of canonical fields ip = 

pl/2 e iS/h ig nQt at aU 

surprising (and rather desirable!). Yet the canoni- 
cal formalism applies equally well to both classical and quantum ensembles. 
Imposing a restriction to linear constraints, for quantum ensembles only, is 
thus rather difficult to justify on physical grounds, and the full set of possible 
canonical constraints therefore deserves serious consideration. 

4 Superselection via canonical constraints 

It has been shown that the set of canonical constraints for configuration- 
space ensembles includes those corresponding to nonlinear constraints on the 
wavefunction. Thus, although solutions of the Schrodinger equation Eq. © 
are closed under arbitrary superpositions, the presence of such constraints in 
general acts to rule out a number of such superpositions, leading directly to 
superselection-type rules. 

Several examples are given here, related to superselection of energy, spin 
direction, and 'classicality'. Note that these examples are intended only to 
be indicative of the manner in which nonlinear constraints can lead to state 
superselection rules. Further investigation is needed to determine whether 
well-motivated constraints exist that lead to superselection rules of physical 
interest. 
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4.1 Energy superselection 

Consider first the rather simple canonical constraint 

J:=PV5 = 0. (9) 

This constraint is local, invariant under S — > S + constant, and may be phys- 
ically interpreted as the requirement that the ensemble momentum density 
vanishes everywhere (see also appendix). Note that for quantum ensembles 
it can be re-expressed as lmijj*'Vijj = 0, which clearly cannot be put in the 
linear Dirac form of Eq. (JJJ). 

To investigate this constraint for a classical ensemble of particles, note 
that consistency with the equations of motion requires dJ/dt = |9*1 IT51 ITS]. 
Eqs. (jHJ) and (0) then yield the secondary constraint 

= d{PVS)/dt = {dP/dt)VS + P(d(VS)/dt) 

= ~[V.(PVS)]VS -PV[\VS\ 2 /(2m) +V] 

= -PVV. (10) 

Hence the classical force, — VV, vanishes over the suport of the ensemble. 
Note in particular that if the potential energy has a single minimum, then the 
constraint requires the ensemble to be concentrated solely at this minimum, 
i.e., the ensemble must occupy the classical groundstate. 

In contrast, for a quantum ensemble, Eq. flOJ) requires that VS* vanishes on 
the support of the wavefunction, and hence that S has no spatial dependence 
for P^O. Secondary constraints arising from consistency with the equations 
of motion can be determined similarly to the classical case above (one finds 
that the 'quantum' force, — V(V + Q), must vanish over the ensemble, where 
Q is the quantum potential in Eq. JSJ)). However, it is simpler to directly 
substitute the ansatz S(x,t) = — fit) into the Schrodinger equation Eq. (0) 
and use the continuity equation (J3J), to obtain the secondary constraints 



fpl/2 



-n 



2 



2m 



■V 2 + V 



P 1/2 , dP/dt = 



respectively. Differentiating the first of these with respect to time and apply- 
ing the second implies that f = E = constant, and hence these constraints 
are equivalent to the time-independent Schrodinger equation 

-h 2 

V 2 ^ + Vtb = Eip. (11) 

2m 
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Thus the quantum ensemble is required to be in an energy eigenstate (but 
not necessarily in the quantum groundstate). 

It is seen that in both the classical and quantum cases, the primary con- 
straint in Eq. © leads to the requirement that the ensemble is stationary. 
In the quantum case this immediately yields a state superselection rule: su- 
perpositions of states of different energy are forbidden. Thus this constraint 
provides a very simple example of how canonical constraints can lead to 
superselection-type rules for quantum ensembles. 

Note that for degenerate Hamiltonians the above example is even more 
stringent than usual superselection rules. In particular, if ipi and ^2 are two 
real-valued energy eigenstates corresponding to the same energy eigenvalue 
E, then it follows from Eq. ® that only superpositions of the form 

e<+ (0^1 + H2) , 

are permitted, where a and b are real. Thus the physical states are not only 
restricted to members of a set of orthogonal subspaces of the Hilbert space; 
they are further restricted, up to an arbitrary phase factor, to real super- 
positions of (real-valued) members of these subspaces. The constraint thus 
effectively lowers the degeneracy of the Hamiltonian (it would be interesting 
to investigate the thermodynamic implications). 



4.2 Spin-direction superselection 

As an example involving a discrete configuration space, an ensemble of two- 
level systems will be considered. The configuration space is therefore labelled 
by two values, j = 1,2, and the canonical equations of motion have the form 
(see appendix) 

S--™ (12) 

for some ensemble Hamiltonian H(P\, P2, Si, S2). 

For the particular example of a quantum ensemble of spin-half particles, 
in an external magnetic field B, the ensemble Hamiltonian has the form (see 
appendix) 

H(P 1 ,P 2 ,S 1 ,S 2 ) := n(a.B) 

= fi(P 1 ~P 2 )B z 

S\ — S 2 . Si — s 2 

B x cos B v sin 

h h 



+ 2fiJPiP 
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where /x is the magnetic moment; a denotes the Pauli sigma matrices; and 
B x , By and B z denote the components of B relative to some set of orthogonal 
directions x, y, z. The canonical coordinates are related to the usual Bloch 
sphere representation by Pi = cos 2 9/2, S1 — S2 = 4>, and it is well known that 
evolution corresponds to rotation of the Bloch sphere about an axis parallel 
to the magnetic field B. Such rotations thus correspond to a particular group 
of canonical transformations of the canonical coordinates Pi, P2, Si and S2 ■ 
Consider now the constraint that 

St - S 2 = (13) 

for some choice of x, y, z. This constraint is analogous to Eq. (0), and is sim- 
ilarly invariant under the transformation Sj — > Sj + C (and hence consistent 
with the conservation of probability). It restricts the ensemble to lie on a 
great circle on the Bloch sphere (passing through the z-axis), and thus ge- 
ometrically corresponds to a rather natural geodesic constraint in the Bloch 
representation. Note, however, that it cannot be expressed as a Dirac-type 
constraint on the wavefunction. 

The secondary constraints arising from the requirement of compatibility 
of Eq. (fT3*j) with the equations of motion follow directly from geometrical 
considerations: rotation of the Bloch sphere about B will leave the ensemble 
on some great circle if and only if (i) the axis of rotation is orthogonal to the 
great circle, or (ii) the ensemble lies on the intersection of the great circle 
with the axis of rotation B. 

The constraint thus leads to a state superselection rule for spin direction: 
the ensemble either corresponds to one of the two eigenstates of cx.B, or to an 
equally-weighted-modulus superposition of these eigenstates. Note that this 
superselection rule is consistent with measurement of spin via a Stern-Gerlach 
apparatus, which leaves the spin-direction either parallel or antiparallel to 
the magnetic field of the apparatus. 

4.3 Gaussian superselection 

The ensemble Hamiltonians H c and H q , for classical and quantum ensembles 
of particles respectively, differ by a functional of the position probability 
density P, as indicated in Eq. (J1J). This functional vanishes in the formal 
limit h — > 0, as one might expect. However, in the real world one cannot 
take such a limit - H is a fixed and fundamental constant - and hence any 
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physical description of a classical limit for quantum ensembles must take a 
different approach [TU] . 

Here one very simple criterion for quantum ensembles to behave classi- 
cally is considered: the difference between the quantum and classical ensemble 
Hamiltonians should be as small as possible. It follows from Eq. (j3j) that this 
criterion is equivalent to the constraint 

F := J dx P |V log P\ 2 = minimum. (14) 

Note that this cannot be expressed as a Dirac-type constraint on the wave- 
fuction. 

The quantity F in Eq. (fTf]) is just the (classical) Fisher information of P 
[2*U| . It is strictly positive, and satisfies the isoperimetric inequality 

F > (27ien)e- 2Hx/n , 

where Hx denotes the entropy of the ensemble and n is the dimension of 
the configuration space, with equality holding if and only P is Gaussian 
[20J. It follows immediately, under the additional (rather weak) assumption 
that the entropy of P is finite, that the classicality constraint in Eq. (fPfj) is 
automatically satisfied if P can be chosen to be Gaussian, i.e., of the form 

P(x, t) = (27rer n/2 (det Kt y^ e -W^-t) T K t (x~a t ) (15) 

for some (possibly time-dependent) non-singular matrix K t and mean a t . 
Note that the finite entropy assumption holds, for example, whenever the 
configuration space variances Var(xi), Var(x 2 ), ... are finite. 

Substitution of Eq. ()15j) into Eqs. (0) leads in a straightforward manner 
to the secondary constraints that each of 5* and V are at most quadratic in x. 
Hence, for quadratic Hamiltonians (eg, the harmonic oscillator, a particle in 
a constant magnetic or gravitational field, and free particles), the classicality 
constraint (|14|) requires the wavefunction to be a complex Gaussian in x. 

The classicality constraint thus leads to a rather strong state superselec- 
tion rule for systems with quadratic Hamiltonians (and finite position en- 
tropy): the wavefunction is restricted to a set of non- orthogonal states. The 
connection between such 'classical' wavefunctions and classical trajectories 
has been recently reviewed in |21j . 
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5 Conclusions 



State superselection rules and nonlinear constraints on the wavefunction are 
inextricably linked in a formal sense. For example, any superselection rule 
of the form that physical wavefunctions are restricted to a set of orthogonal 
subspaces, corresponding to an orthogonal set of projections {Ej} (and hence 
to eigenstates of a Hermitian operator A := J2j a jEj), is formally equivalent 
to the nonlinear constraint 

£<^#> 2 = 1 (16) 

3 

on the wavefunction. It follows that if state superselection rules have physical 
content, then so must nonlinear constraints, and vice versa. 

It has been demonstrated, in the context of configuration-space ensembles 
satisfying an action principle, that the consideration of nonlinear constraints 
is in principle unavoidable. Moreover, in this context both quantum and 
classical constraints can be discussed on an equal footing, generalising the 
standard approach based on mapping classical phase space constraints to 
linear operator constraints. 

The examples of constraints leading to superselection of energy, spin di- 
rection, and 'classicality' in section 4 show that quite simple nonlinear con- 
straints on the wavefunction can lead to interesting state superselection rules. 
However, further investigation is needed to determine whether well-motivated 
constraints exist (i.e., not of the ad hoc form of Eq. (Jrfijl above), that lead to 
superselection rules for quantities such as total charge and spin discussed in 
the Introduction. It would also be of interest to investigate what conditions 
need to be imposed, if any, to ensure the consistency of such superselection 
rules with physical measurements. 

The canonical formulation of the equations of motion emphasises the fact 
that an action principle for quantum ensembles requires two fields. These 
are usually taken as if> and if>* j7j, but of course can also be taken as P and 
S. Interpretations of quantum mechanics which neglect this fact (i.e., which 
neglect to provide a physical interpretation for both if) and if}*, or both P 
and S), correspond to according the quantum equations of motion primary 
significance over the action principle from which they follow. This is opposite 
to the usual point of view taken in classical dynamics [TBI . 

It is perhaps worth pointing out that, since the wavefunction can be re- 
covered from tomographic measurements (up to a phase factor), the fields 
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P and S can similarly be recovered (the latter up to an additive constant). 
Indeed, weak measurements of momentum on an ensemble of particles, post- 
selected by strong measurements of position, allow one to recover both P{x) 
and VS 1 for both the classical and the quantum case [22] (where the latter 
can be integrated to give S up to a constant). 

Finally, note that the canonical formulation of the equations of motion 
only applies to pure states in the quantum case. It would be of interest, for 
thermodynamic and other purposes, to consider properties of statistical mix- 
tures of configuration-space ensembles (corresponding to density operators 
in the quantum case). 
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A Appendix: Configuration-space ensembles 

To introduce configuration-space ensembles at a fundamental level, suppose 
that (i) the configuration of a physical system can only be specified impre- 
cisely, requiring it to be described by an ensemble on the configuration space, 
and (ii) there is a action principle for the motion of the ensemble. 

The first assumption implies that the configuration of the physical system 
is formally described by a probability density P on configuration space. The 
second assumption implies the existence of some conjugate function S on 
configuration space, and a 'ensemble' Hamiltonian H[P,S,t], such that the 
equations of motion are generated by Hamilton's action principle 5a = 0, 
with [E] 

Jdt H + J P^j- = Jdt[H+ (dS/dt)] . (17) 



a 



The inner integral in the first equality above is over the configuration space, 
and is replaced by summation over any discrete parts of configuration space. 
The physical interpretation of S will be discussed further below. For contin- 
uous configuration spaces it is typically an energy-momentum potential. 

It is perhaps simplest to first consider discrete configuration spaces (corre- 
sponding, eg, to the description of classical dice throws and quantum spins). 
The probability of configuration j is then Pj, with conjugate quantity Sj, 
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and the equations of motion follow from Eq. (J 17)) as 

Pj = OH/dS j} Sj = -dH/dPj. (18) 

Only Hamiltonians which conserve probability are of interest, and hence, to 
first order in e, it is required that 

arr 

= eZPj = eZ%- = H(P,S + e)-H(P,S). 

3 3 ° b 3 

Thus, H is invariant under the addition of a constant to the components of 
S. 

It follows for discrete configuration spaces that the ensemble Hamiltonian 
must be of the form 

H = F(P,M) } (19) 

where M denotes the antisymmetric matrix with components Mjk = Sj — Sk- 
The equation of motion for Pj therefore reduces to the familiar rate equation 
form 23 

^=?te-S=? (r * a - r ^ )> (20) 

with transition rates 

T ^ l m (21) 

(with Tjk := for P*. = 0). Noting that probability must remain positive 
under evolution implies further that 



OH 



dSj 



P j= k&j) 



> 0. (22) 

Pi=0 



For constant transition rates this reduces to 7}^ > for j ^ k. 

For a continuous configuration space, indexed by some (possibly multi- 
dimensional) parameter x, the ensemble Hamiltonian becomes a functional, 
H[P, S], and the equations of motion follow from Eq. (fT7|) as 

dP__5H_ dS_ _SH_ 

'dt'Js' ~dt~~~5P' 
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Here SF/Sf denotes the functional derivative of functional F[f] with respect 
to function /, where for the special case that F[f] = J dx R(f, V/, x) one 
has PH 

SF/Sf = dR/df -V.[dR/d{Vf)). (23) 

The conservation and positivity of probability places restrictions on the 
ensemble Hamiltonian analogous to the discrete case above. In particular, 
conservation of probability implies that H is invariant under S(x) — > S(x) + 
constant, and so can only depend on VS and higher derivatives. The equation 
of motion for P can therefore always be expressed in the form of a continuity 
equation, 

dP/dt + V.J = 0, 

for some current J (one formally has J = SH/S(VS)). 

Examples of ensemble Hamiltonians for classical and quantum particles 
have been given in Eqs. (J2J) and (@J) respectively (with current J = PVS/m in 
both cases). The form of the classical Hamiltonian H c is perhaps the simplest 
possible scalar functional of VS 1 , and is well known (dating back to the 19th 
century) from the theory of ideal fluids [21]. The form of the quantum 
Hamiltonian H q may be derived from the corresponding classical form via 
an 'exact uncertainty' principle, both for particles jT5j and for bosonic fields 

m- 

It is straightforward to show that all quantum systems may be formulated 
in terms of configuration-space ensembles, where the 'configuration space' 
may be chosen as corresponding to the outcome space of any complete mea- 
surement on the system. In particular, let H be the Hamiltonian operator 
for a quantum system described by state \ip), and let {\a)} denote any basis 
set satisfying the completeness property 

J da | a) (a\ = 1 

(with integration replaced by summation over any discrete ranges of a) . Such 
a basis set corresponds to the outcomes associated with the measurement of 
some observable A (eg, position, momentum, optical phase, energy, spin, 
etc.), where |(a|?/>)| 2 is the probability density associated with measurement 
outcome a. Note it is not required that the basis set be orthonormal. 

Now define P(a) and S(a) by the polar decomposition (a\vp) = p 1 / 2 e tS / h , 
and the associated ensemble Hamiltonian Ha by 

H A [P,S] = (ip\H\iP) = J dada'[P(a)P(a')] 1/2 (a\H\a')e i[s{a ' ) - s{a)]/n . (24) 
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It follows that the action in Eq. (|17|) can be rewritten as 



a — dt 




am _ m -?m 

and hence (by variation with respect to (a\ip) and its conjugate) that the 
equations of motion for P and S are identical to the Schrodinger equation 
(and its conjugate) in the a-representation. 

The natural form of the transition rates for a discrete quantum observable 
A follows from Eqs. (|2TJ) and (J2IJ) as 

T jk = h-\P /P k ) l l 2 \m {{ aj \H\a k )e-^- s ^ h } . (25) 

The positivity condition (}2*2*|) automatically follows (with equality) from the 
Hermiticity of H. An alternative 'geometric' approach to the configuration- 
space ensemble formulation, for the case of discrete quantum ensembles, has 
been given recently by Mehrafarin jT^j. 

It is of interest to give an explicit example of the quantum equations of 
motion for the case of a non- orthogonalhasis set (see also j2H])- I n particular, 
consider a single-mode bosonic field with Hamiltonian operator H = TiojN . 
The phase observable <£> then corresponds to the basis set {\<fi)} with number- 
state expansion l2H] 



oo 

|0) ;= (2vr)- 1 / 2 ^e^|n) 

n=0 



It follows that {(p\N\ijj) = i(d/d(p){<j)\if)), and hence that the corresponding 
ensemble Hamiltonian is given by 

HJP,S] = hu}(if>\N\iJ)) = -u / d<j)P—. 

Jo o<p 

The associated equations of motion follow as 

dP/dt = u (dP/d(p) , dS/dt = u (dS/d<f)) . 

These may be solved directly, or by noting that H$ is proportional to the 
ensemble momentum associated with phase translations (cf. Eq. (J28)) below), 
to give P((j>, t) = P o (0 + ut) and 5(0, t) = S (<p + ut). 



16 



Finally, to establish the physical interpretation of the conjugate quantity 
S, recall from the above discussion that the conservation of probability im- 
plies invariance of H under the transformation S — > S + C. The equations of 
motion are also invariant under this transformation. Hence only derivatives 
and relative changes in S are expected to have direct physical significance. 

For example, the form of the ensemble action in Eq. (J 17)) suggests that 
dS/dt is a local energy associated with the ensemble. This can be made 
more precise, both for continuous and discrete configuration spaces, in the 
case that the ensemble Hamiltonian is homogeneous with respect to P, i.e., 

H[XP,S]= X r H[P,S] (26) 

for all A > and some constant r. In particular, differentiating the homo- 
geneity condition with respect to A and putting A = 1 yields the identity 

H = r~ l [ dxP^- = -r- 1 [ dx P^- = -r^ldS/dt) (27) 
J oP J ot 

(integration is replaced by summation for discrete configuration spaces). 
Thus, identifying H with the ensemble energy (which is certainly justified 
for time-independent Hamiltonians) , it follows that —r~ 1 P{dS/dt) may be 
interpreted as a local energy density. This is further supported by the invari- 
ance of the equations of motion under H — > H + e(J dxP) r , S — > S — ert. 
Note that the homogeneity property (}2T)j) in fact holds for all cases considered 
in this paper, with r = 1. 

Further, for continuous configuration spaces, consider an infinitesimal 
translation in configuration space. The change in the ensemble Hamiltonian 
is then, to first-order, 

H[P(x + 6x),S(x + 6x)]-H[P(x),S(x)] = Sx. J dx VP + || Vs) 

= Sx.i- I dxPVS, 



dt 

where the last equality follows using the equations of motion and integration 
by parts. It follows that the linear momentum of the ensemble may be 
identified as [in] 

n = J dxPVS = (VS), (28) 

and hence that PV5* may be interpreted as a local momentum density asso- 
ciated with the ensemble. 
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